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Abstract. A linear theory of whistler wave is developed wihtin the paradigm of 
a two dimensional incompressible electron magnetohydrodynamics model. Exact 
analytic wave solutions are obtained for a small amplitude whistler wave that ex- 
hibit magnetic field topological structures consistent with the observations and our 
simulations in linear regime. In agreement with experiment, we find that the par- 
allel group velocity of the wave is large compared to its perpendicular counterpart. 
Numerical simulations of collisional interactions demonstrate that the wave mag- 
netic field either coalesces or repels depending upon the polarity of the associated 
current. In the nonlinear regime, our simulations demonstrate that the evolution of 
wave magnetic field is governed essentially by the nonlinear Hall force. 



1. Introduction 

Whistler waves are widely observed in many space and laboratory plasma phe- 
nomena. For instance, they are believed to be generated in the Earth's ionospheric 
region by lightning discharges and proceed in the direction of Earth's dipole mag- 
netic field [1]. They have also been recently detected in the Earth's radiation belt 
by the STEREO S/ WAVES instrument [2]. Moreover, there are observations of 
Venus' ionosphere that reveal strong, circularly polarized, electromagnetic waves 
with frequencies near 100Hz. The waves appear as bursts of radiation lasting 0.25 
to 0.5s, and have the expected properties of whistler-mode signals generated by 
lightning discharges in Venus' clouds [3]. These waves are also reported near the 
Earth's magnetopause [4] and the Cluster spacecraft encountered them during the 
process of magnetic reconnection in the Earth's magnetotail region [5]. Upstream 
of collisionless shock, whistler waves are found to play a crucial role in heating the 
plasma ions [6] . Their excitation and propagation are not only limited to the Earth's 
nearby ionosphere, but they are also found to be excited near the ionosphere of other 
planets such as in the radiation belts of Jupiter and Saturn [7] . Whistler waves are 
believed to be a promising candidate in transporting fields and currents in plasma 
opening switch (POS) devices [8], which operate on fast electron time scales. Simi- 
larly, these waves have been known to drive the phenomenon of magnetic field line 
reconnection [9] in astrophysical plasmas [10]. Whistler waves have also been inves- 
tigated in several laboratory experiments [11, 12, 13, 14, 15, 16], where they have 
been found to exhibit a variety of interesting features such as anisotropic prop- 
agation of phase front, strong dispersion characteristics, interaction with plasma 
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particles etc. A few experimental features have also been confirmed by recent three 
dimensional simulations [17], where it has been reported that the polarity and the 
amplitude of the toroidal magnetic field, in agreement with the laboratory experi- 
ments, determine the propagation direction and speed of the whistlers. These are 
only a few examples amongst a large body of work devoted to the study of whistler 
waves. Despite the large amount of efforts gone into understanding the existence 
and propagation of whistlers, their linear dynamics is still debated, specially in the 
context of complex nonlinear processes. For example, the role of whistler waves 
in high frequency turbulence and anisotropic spectral cascades has been debated 
recently [18, 19, 20, 21, 22]. 

Motivated by these issues, we in this paper develop a fully self-consistent linear 
analytic theory of whistler wave propagation. The electromagnetic whistler waves in 
a plasma are excited essentially by means of collective electron oscillations in pres- 
ence of an external or self-consistent large scale magnetic field. Their characteristics 
frequency thus lies between ion- and electron-gyrofrequencies i.e., uj C i <C u) <C w ce , 
where ui C i and u ce are ion- and electron-gyrofrequencies respectively. In such a fre- 
quency regime, usual magnetohydrodynamic theory becomes invalid and the plasma 
processes need to be described by electron magnetohydrodynamics (EMHD) theory 
[23]. The linear properties of EMHD are mainly governed by whistler waves, which 
constitute fundamental oscillatory modes of this model. 

In this paper, we develop a linear analytical theory and simulation to investigate 
the linear propagation of a whistler wave packet and demonstrate that such waves 
lead to an exact solution of EMHD model. For the sake of simplicity, we restrict 
ourselves to a two dimensional plane where variations in the wave magnetic field arc 
confined to the x and y directions only. The analytic results of our theory show a 
qualitative/quantitative agreement with the experimental observations and simula- 
tions. Such a comparison further enables us explain several experimentally observed 
features of the propagating whistler wave packets. We also carry out an investiga- 
tion of collision of two propagating whistler wave packets to explore their mutual 
interaction in the context of ionospheric processes. The remainder of this paper is 
organized as follows. In Sec 2, governing equations of whistler waves and dispersive 
properties are described. In Sec 3, we describe the evolution of wave amplitude 
based on Fourier analytic method and numerical simulations. Section 4 describes 
the propagation studies of whistlers that show a remarkable agreement with the 
simulation in the linear regime. The linear simulations are discussed in Sec 5. Sec- 
tion 6 contains results of mutual interaction of two whistler modes, whereas Sec 7 
describes nonlinear evolution of the whistlers. Finally conclusions are discussed in 
Sec 8. 

2. Dynamical Equations of Whistler Waves 

Whistler wave dynamics is governed essentially by the equations of EMHD. The 
EMHD phenomena occur typically on rapid electron time scales [23], while ions 
do not participate in the dynamics. Thus, the basic frequency scales involved are 
ujd <C u> -C oJce (where u c j, lu C€ are respectively ions and electrons gyro-frequencies, 
and lo is the characteristic frequency) and the length scales are c/uo pi < I < c/uo pei 
where w p j, u pe are the plasma ion and electron frequencies. Currents carried by the 
electron fluid are important. Since stationary ions merely provide a neutralizing 
background to a quasi-neutral EMHD plasma, their momentum can be ignored in 
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the high frequency regime. The momentum equation of electron is 

m e n ( + V e • VV e J = -enE - — V e x B - VP - ^m e nV e , (2.1) 

1 f)A 

E^-V^--, (2.2) 
VxB = ^J + if, (2.3) 

^+V-(nV e )=0. (2.4) 

The remaining equations are B = V x A, J = — enV e , V • B = 0. Here m ei n, V e 
are the electron mass, density and fluid velocity respectively. E, B respectively rep- 
resent electric and magnetic fields and <j>, A are electrostatic and electromagnetic 
potentials. The remaining variables and constants are, the pressure P, the colli- 
sional dissipation /n, the current due to electrons flow J, and the velocity of light c. 
The displacement current in Ampere's law Eq. (2.3) is ignored, and the density is 
considered as constant throughout the analysis. The electron continuity equation 
can therefore be represented by a divergence-less electron fluid velocity V • V e = 0. 
The electron fluid velocity can then be associated with the rotational magnetic field 
through 

V e = — V X B. 

Anne 

On taking the curl of Eq. (2.1) and, after slight rearrangement of the terms, we 
obtain a generalized electron momentum equation in the following form. 

^ - V e X (V x P) + VC = -M™ e V e (2.5) 

where 

eA 1 P 

P = m e V e - — and £ = -m e V e ■ V e H e<p. 

c 2 n 

The curl of Eq. (2.5) leads to a three-dimensional equation of EMHD describing 

the evolution of the whistler wave magnetic field, 

-|(B - d 2 V 2 B) + V e ■ V(B - dlV 2 B) - (B - ^ V 2 B) • VV e = ^d 2 e W 2 B. (2.6) 

where d e ~ c/uo pei electron skin depth, is an intrinsic length-scale in EMHD. The 
three-dimensional equations of EMHD can be transformed into two dimensions by 
regarding variation in the z-direction as ignorable i.e. d/dz = 0, and separating the 
total magnetic field B into two scalar variables, such that B — z x V-0 + bz. Here 
tp and b respectively present perpendicular and parallel components of the wave 
magnetic field. The corresponding equations of these components can be written in 
a normalized form as follows, 

-(V - d 2 e V 2 yj) + zxVb-V(yj- ^V 2 V) - B —b = 0, (2.7) 

— (b - dlM 2 b) - d 2 J x V6 • VV 2 6 + z x V^ • VV 2 ^ + ^o^V 2 ^ - 0. (2.8) 

The length and time scales are normalized respectively by d e and oj ce , whereas 
magnetic field is normalized by a typical mean B . The linearization of Eqs. (2.7) 



4 



Dastgeer Shaikh 



& (2.8) about a constant magnetic field Bo yields the dispersion relation for the 
whistlers, the normal mode of oscillation in the EMHD frequency regime, and is 
given by 

(ig ky k 



ca l + dlk 2 ' 

where uj Co = eB /mc and fc 2 = fc 2 + ky. From the set of the EMHD Eqs. (2.7) 
& (2.8), there exists an intrinsic length scale d e , the electron inertial skin depth, 
which divides the entire spectrum into two regions; namely short scale (kd e > 1) 
and long scale (kd e < 1) regimes. In the regime kd e < 1, the linear frequency of 
whistlers is LOk ~ k y k and the waves are dispersive. Conversely, dispersion is weak 
in the other regime kd e > 1 since tot ~ k y /k and hence the whistler wave packets 
interact more like the eddies of hydrodynamical fluids. 



3. Evolution of Whistler Wave Amplitude 

The evolution of the amplitude associated with a linear whistler wave can be studied 
from Eqs. (2.7) & (2.8) after dropping the nonlinear terms, as described in the 
following. 

V 2 V)-Bo^=MV 2 V (3.1) 

^(6-V 2 6) + B |-VV-MV 2 6. (3.2) 

The rhs of the above equations now contains a damping term. We adopt following 
definition to Fourier transform the evolution variables of Eqs. (3.1) & (3.2), 

G(x,y,t) = — / G(k x ,k y ,t)exp(ik x x + ik y y)dk x dk y (3.3) 

z ™ J — oo J — oo 

and the corresponding inverse Fourier transform is 

G(k x ,k y ,t) = — / / G(x,y,t)exp(-ik x x - ik y y)dxdy. (3.4) 

^ J — oo J — oo 

We then write Eqs. (3.1) & (3.2) as, 

/ -fj,k 2 ikyBpk 2 \ 

/ T+F l+fc^ \ 




iky Bp -fik 2 

\ 1+F T+F / 




(3.5) 



where bk ~ b(k x ,k y ,t) and Vfc ~ il>{k x , k y ,t). The solutions of Eq. (3.5) can be 
written as; 

h(t) = \ MO) + fc^fe(O)] exp(A+t) + 1 [& fc (0) - fcVfe(O)] exp(A_t) (3.6) 

Mt) = ^ MO) + ki> k (0)} exp(A+i) - ±- [6 fe (0) - fc^(O)] exp(A_t) (3.7) 

where 6^(0) and ipk{0) denote the Fourier transformed initial conditions of the 
variables b(x, y, t = 0) and ip(x, y, t = 0) respectively, and A is an eigen value of the 
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squared matrix of Eq. (3.5) which is given as below; 

\ + ik y kB ( o q\ 

A± = TTp ± T+W (3 ' 8) 

Here also, ± sign corresponds to the wave moving in the forward and backward 
direction. To seek an analytic solution of Eqs. (3.6) & (3.7), we initialize the wave 
magnetic field perturbations by means of a current carrying antenna similar to that 
of Stenzil's et al [15] experiment such that there exists only in plane component of 
the wave magnetic field, whereas the component along the ambient magnetic field 
is zero initially. Accordingly, we choose the following kind of initial perturbation to 
study the behaviour of wave amplitude 

/ y 2 

ip(x, y, t = 0) = cos(aa;) exp — — ^ 

V a v 

b(x,y,t = 0) = 
Fourier transformation yields, 

^fc(O) - ^V^[S(k x +a)+ S(k x - a)] exp (-^^j (3.9) 

Mo) - o 

We would like to concentrate on the case in which a wave packet is propagating in 
the forward direction only, and hence choose the eigen value that corresponds to A+ 
in the Eq. (3.8). Under such conditions, the poloidal component of wave magnetic 
field i.e., Eq. (3.7) assumes the form as; 

i>k(t) = ^V'fc(0)exp(A + t) 
substituting Eq. (3.9) into above eqn and then inverse Fourier transformation gives 

V>(x, y, t) = ^ a(J y^ f + [S(k x +a) + S(k x - a)} e^ V {ik x x)Q(k x )dk x (3.10) 

^ J-oo 

where, 

r+oo / 0.2^,2 \ 

C{k x ) = J exp f ^ j cxp(iB kkyt)cxp(ik y y)dky 7 (3.11) 

The expression of ip contains another integral given by Eq. (3.11). We therefore 
first solve Eq. (3.11) and substitute the solution into Eq. (3.10). In order to do so, 
we put Eq. (3.11) into the following form 

/+oo 
F(k) exp [i X (k)t] dk y (3.12) 
-oo 

where k 2 = {k x + k 2 y ) and 

F(k) = exp 



X(k) = uj{k) + k y - 
u)(k) = Bokyk 
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We employ the method of stationary phases [24, 25] to analyze Eq. (3.12). This 
method indicates that in the time asymptotic limit (i.e., at large times), the dom- 
inant contribution in the integral will essentially come from the maximum of the 
exponential part. This further suggests that the function x(fc) in Eq. (3.12), must 
take on a real and positive maximum value. Away from this maximum, the in- 
tegrand will have negligibly small contribution. The condition of maximization of 
x(k) further implies that 

A x(fc) = o. (3.13) 

The vanishing of first derivative reveals that we do have stationary points, which 
may be either maximum or minimum. Equation (3.13) gives a bi-quadratic algebraic 
equation in k y , whose roots are basically functions of k x . 

\ 1/2 



k <k)-±l J?— - M ± y I y2 



8B 2 t 2 2 2B t]j 1652*2 2 

The above equation admits four roots. We thus have four stationary points corre- 
sponding to these four roots. We take into account only that point, which gives a 
maximum contribution in the Eq. (3.12). We denote this point as fco = ^^(kx). 
Taylor expanding the function x(fc) around this stationary point (fco), we may write 

X (k) = x(fco) + (k y - k ) X '(ko) +hk y - k f X "(k ). 

V v ' 2! 



Note the second term in the expansion is zero because of the condition of extremum 
imposed by Eq. (3.13). Substituting x(fc) in Eq. (3.12), we get 



f + OO 

C(k x ) = F(fc ) exp [ix(k )t] I exp 



i^(ky-k ) 2 x"(k )t 



•Mko)t] J 

Equation (3.10) with the help of above expression give, 

i>(x, V, t) = — ^ / [S(k x + a) + S(k x - a)] exp(ik x x) 



dky . 



/ 27T 

xF(fco) exp [i X (ko)t] J tx „( ko) dk x (3-14) 



where we have used exp [i^(k y — k n ) 2 X " {ko)t\ dk y = \J tx n^ ko y The solution 
can be written as 



1>{x, y, t) = C os(ax)F(k a )ex V [i X (k )t] J (3.15) 

4tt y tx"{k ) 

Equation (3.15) can be used to predict the behaviour of wave amplitude at a given 
location in space. The various quantities in this equation are as follows; 



fc 



2 / 9 

y y y 



1/2 



8B%t 2 2B t V l6B 2 t 2 2 2 
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Figure 1. The behaviour of linear wave amplitude. Circle shows analytical prediction of 
wave amplitude, which overlaps on the amplitude obtained from numerical simulation of 
Eqs. (2.7) & (2.8) in the linear regime. Clearly, the analytic amplitude and that obtained 
from numerical simulation show an excellent agreement in a long time evolution. 



F(k ) = exp - 



a y K 



X(fco) = B k (k 2 + a 2 ) 1 ' 2 + fc 



X"(M = B 



3fc 



VW+^ (fc 2 + « 2 ) 3/2 



We next numerically integrate Eqs. (2.7) & (2.8) with the help of pseudospectral 
code to determine the amplitude of the propagating whistler wave in the linear 
regime. This amplitude will then be compared with the one predicted by Eq. (3.15). 
The numerical code is based on the discrete Fourier representation in fc-space in 
the two directions viz, x and y. The boundary conditions in the two directions are 
periodic in nature. The linear part of the equations is exactly integrated in Fourier 
space. On the other hand, the evaluation of nonlinear part is carried out in real space 
and then Fourier transform it to fc-space. We use Fast Fourier Transform (FFT) 
routines to go back and forth in the real and fc-space at each time integration step. 
The time advancement is done with the help of a second order leap-frog predictor 
corrector scheme. The numerical accuracy of the results is checked continuously by 
monitoring the energy conservation laws of EMHD equations at each time interval. 
We evolve 64 2 and/or 128 2 Fourier modes in the two dimensional box of size 207T x 
207T and measure the amplitude of wave at a spatial location determined by the 
computational box co-ordinates (x = 9.5n,y = 9.5n). 

A comparison of the amplitude obtained from linear simulation and stationary 
phase method is shown in Fig. (1). It can be seen from the figure that the two results 
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are in quantitative agreement with each other. Such a close agreement between 
analytical and numerical results thus clearly demonstrates the validity of the our 
simulation code in a linear regime. 



4. Propagation Theory of Whistler Wave - Exact Wave Solutions 

We now investigate the case of a propagating whistler wave packet, which is excited 
by a current-carrying linear wire antenna. Such a study is motivated primarily by 
the work described in [1-10] and specifically by the experimental observation of 
Ref. [15]. In these experiments [15], it has been observed that the linear whistler 
waves can be excited in EMHD plasma, when a current-carrying linear antenna 
is kept across an external magnetic field. The propagating wave packets assume 
conical shape whose parallel group velocity is more than the perpendicular one. 
Furthermore, the wave propagation has not been observed in their experiment, 
when the exciter was placed along the direction of an external magnetic field. In 
order to understand the propagation characteristics of whistler waves, we carry out 
both analytical as well as numerical calculations in the linear regime. We first study 
the evolution of wave pattern using analytical treatment. 
The initial conditions are as follows; 



b(x,y,0) = 



(4.1) 



where V^o and a are amplitude and width of the distribution of the initial wave 
packet respectively. Equation (4.1) represents constant circular (if>) contours of wave 
magnetic field in the plane of variation (i.e., in the xy-planc). Such perturbations 
in the xy-pl&ne can arise as a result of a current-carrying wire antenna, which is 
kept parallel to the z-direction. Hence there is no perturbed magnetic field along 
this direction. The initial condition for the axial component of the wave magnetic 
field is chosen to be zero, b(x, y, t = 0) = 0. 

4-1. Analytical Treatment : 

We Fourier transform Eq. (4.1) using the definition of Eq. (3.3) and substitute it 
in Eq. (3.7) to get, 

, , . * o" 2 ( o 2 k 2 \ , , ... 

Vk{t) = — - — cxp — {exp(A + t) + cxp(A_i)} 



We inverse Fourier transform this equation and solve in the limit of kd e <C 1 and 
/j, = 0. Under these limits, the eigen values are A± = ±iB k y k, and the above eqn 
becomes tp(x, y, t) = I + + 7_, where 

^ a 2 r+oo I- + 00 / 0.2^,2 \ 

I± = — — / I dk x dk y exp(ik x x + ik y y) exp ( — ] cxp(zkiBok y kt) 



oo »/ — oo 



We first solve I + . Using the transformations k x = kcos9 and k y = fcsin#, which 
give dk x dk y = k dk d6. The angle 9 varies from to 27r, while k varies from to 
oo. The integral I + then takes the form, 

* cr 2 f +oc f 27T ( <r 2 k 2 \ 
I + = — — J j kdkdO cxp ^ — J exp [i(ky + B a k 2 t) sin#] 
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Figure 2. Propagation of whistler waves. The wave fronts have conical shape. Figure shows 
constant i/>-contours at various time. The vertical axis represents j/-axis, along which Bo 
is oriented. Horizontal axis indicates a>axis. Similar wave field topology has been reported 
inRef[15]. 

x exp (ikx cos 9) 

The exponential term containing sine in the above expression can be expanded 
using Bessel identity 26, exp(i£sin0) = J2n=-oo Jn(0 exp(m0). We thus get 



1+ 



8tt 



2 r+oo 
JO JO 



kdkdQ exp 



Jn(ky + B Q k 2 t) 

— — oo 

x exp [i(n9 + kx cos 9)] 



We then carry out the 9- integral with the help of the identity [26], exp[i(n9+ 
(3 cos 9)]d9 — J„(/3), and get the equation, 



2 r+oo 

/ kdk exp 
Jo 



^ J n (ky + B Q k 2 t) exp 



inir\ 

— J Jn(kx) 



Using one of the Bessel identities [26], we then evaluate the summation of above 
equation. This gives us 



U = 



*oo- 2 



2 r+°° 



kdk exp 



a 2 k 2 



o (^(xkf + (yk + B n kH) 2 ) 



We carry out the similar analysis to integrate /_ and get the total ip as follows, 



[j (n+) + j (n_ 



(4.2) 



Similarly, we solve for Eq. (3.6) and get the expression for the axial component as 
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Figure 3. Constant 6-contours of magnetic field of propagating whistler waves. Note that 
the b component is zero initially, and is generated later by means of linear mode coupling 
process. 



follows, 



b(x,y,t) = 



* cr 2 



dfcexp 



o 2 k 2 



[j (n+)- j (n_)] 



(4.3) 



where H± = (xk) 2 + (yk ± B^kH) 2 . Equations (4.2) & (4.3) form an exact whistler 
wave solution of the EMHD model representing respectively the evolution of poloidal 
and axial components of wave magnetic field. The substitution of t = in Eqs. (4.2) 
& (4.3) retrieves the initial condition i.e., Eq. (4.1). The numerical integration of 
Eqs. (4.2) & (4.3) at various time are displayed in Figs. (2) & (3). It is clear from 
the figures that the propagating wave fronts are conical in shape. It can be noticed 
from Fig. (3) that the axial component, though absent initially (i.e., at t = 0), 
acquires a finite amplitude during propagation. 

The axial component of wave field is generated primarily due to linear coupling 
in the equations representing the evolution of b and tp (i.e., Eqs. (3.2) & (3.1)). The 
topology of the magnetic field in the propagating wave packet is consistent with 
Ref [15]. 



5. Simulation of Linear Wave Propagation 

We now carry out numerical simulation of low amplitude whistler waves. The initial 
condition is the same as that used earlier in the analytical treatment. However, 
unlike the linearized approximation (Eqs. (3.2) & (3.1)) used earlier to obtain the 
analytical expression for the evolved fields, here we make use of the exact Eqs. (2.7) 
& (2.8) for simulation. At small amplitude we expect the linear approximation to 
be valid, and hence an agreement of the simulation results with those presented in 
the earlier section. 
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Figure 4. Diffusion of wave magnetic field, when the exciter is parallel to Bo- (a) repre- 
sents constant contours of perturbed magnetic field at t = and (b) shows the diffusion 
of initial field at t = 20 unit of time. 



In the absence of Bq, the fields simply diffuse as shown in Fig. (4). This config- 
uration can be viewed as identical to that of the Ref [15] wherein the linear wire 
antenna is kept parallel to the external magnetic field. The field lines merely dif- 
fuse in such a scenario. The other case of linear wire antenna being placed across 
an external magnetic field can be simulated by choosing Bo as finite. In this case 
we observe propagating whistler wave packets with asymmetric orientation in the 
rry-plane. The phase fronts have conical shape as reported both in our approxi- 
mate linear analytical treatment of the problem as also in experiments. This can be 
understood by realising that the fronts move faster along y (direction of external 
magnetic field) than across. This also shows clearly that the group velocity along 
B is larger than the perpendicular group velocity. A quantitative estimation of 




the group velocities of the propagating whistler wave packet in the two directions 
i.e., parallel and perpendicular to B can be carried out as follows. The two group 
velocities can be computed from the dispersion relation. The perpendicular group 
velocity (Vg ± ) can be given as 

_ dw __ B a kyk x 2kk x k y B 

G± ~ dk x ~ k(i + k 2 ) (i + fc 2 ) 2 ' 1 ' ' 
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The typical wave numbers k x and k y in the above equation can be estimated as a~ x 
and a-y 1 respectively, which represent the width of initial Gaussian distribution of 
the wave packet. Replacing k x and k y by er" 1 and a' 1 respectively in Eq. (5.1), we 
get 

_2Bo 

B /cr x a y (r x <r t 



* v^FWI 1 + (* + *)} {' + (* + *)}' 

Similarly, the parallel group velocity gives, 



(5.2) 



Vg„ 



Bo 



v / S r ^W( 1+ ^ + ^) 



1 Q I ~T~ 



f 7 : 



(5.3) 



We had chosen in the simulation, a% = = 20.0 and Bq = 1.0, which yields 
Vg ± = 0.14 and Vq^ = 0.41, thereby clearly indicating that Vc n > Vq ± - It is 
because of the larger parallel group speed, the wave magnetic field structures are 
conical and they propagate along the direction of the ambient magnetic field. 





Figure 6. Interaction of two whistler wave packets excited by parallel currents in the two 
antennae show an attraction. Figure (a) shows the wave fields at t = and figure (b) 
represent the field at t = 20. 



6. Collisional Interaction of Two Whistlers 

We now study the collisional interaction between whistler wave packets. Such inter- 
actions can possibly exist in the Earth's ionosphere where the oppositely traveling 
waves along the Earth's dipolc magnetic field collide. In a laboratory plasma, such 
an interaction can be studied by exciting the wave field from the two linear wire 
antennae separated by a finite distance. The currents carried by the oppositely 
traveling waves are either parallel or anti-parallel to each other. Our objective here 
is to investigate as what kind of resultant state would be borne out from such in- 
teractions. When initial excitation is due to the antennae current flowing in the 
same direction (i.e., parallel currents in the antennae), they give rise to wave fields 
of similar signs in the plasma. This is shown in Fig. (6) (a). The evolution of the 
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two wave packets then leads to the formation of a magnetic null line. The magnetic 
pressure being weak at the central region, there is a tendency of attraction between 
the magnetic configuration caused by the two currents. Along with this attraction 
as the whistler waves get excited and propagate, the null line at the center alters 
and forms an O-point (see Fig. (6)(b)). The attraction of wave magnetic fields can 
also be understood as follows. The wave magnetic field induces parallel currents in 
the plasma, which then attract each other. This can be contrasted with Fig. (7) 
where the initial excitation are due to the anti-parallel currents in the two wire 
antennae. It is observed that there is no such magnetic null line for this case, albeit 
here the magnetic field maximizes in the centeral region, causing repulsion. It can 
be clearly seen that in this case there is no formation of O-point at the center. This 
is shown in Fig. (7)(a,b). In this case, the wave fields induces anti-parallel currents, 
which repel each other. 

30 F — — 
20 
10 


-10 
-20 

-30^ — — 
-20 

Figure 7. Two whistler wave packets launched by anti-parallel currents in the two wire 
antennae repel each other. Figures (a,b) show the wave fields respectively at t = and 
t = 20. 




7. Nonlinear Whistler Waves 

In the previous section we have presented detailed studies on the propagation and 
collision of whistler waves, when the amplitude of initial perturbation was smaller 
than the external magnetic field (Bo). We now carry out our investigation, when 
the amplitude of the perturbed field is comparable or greater than the ambient 
magnetic field (Bo). This is the regime where nonlinear interactions may be of great 
importance. Nonlinear excitation of whistler in laboratory plasmas has n't been yet 
understood fully [16]. We therefore undertake this study to understand nonlinear 
evolution of whistlers. Our choice of initial field perturbation is identical to the one 
discussed in the previous sections, which was also consistent with the excitation 
mechanism adopted in the experiments [15, 16]. In this section we choose to study 
the evolution of the similar initial perturbation, however, now the perturbed field 
has an amplitude which is comparable or larger than the ambient field -Bo- Figure 
(8) shows the initial and the evolved field configurations. Clearly the magnetic field 
topology is distorted and distinct from the linear results presented earlier (see Fig. 
(2), Fig. (5)(a)). 

The distortion (or modification) in the propagation characteristics of whistler 
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Figure 8. Large amplitude whistler wave. The wave field is distorted in —x direction 
due to experiencing a net J x B force. 



waves can be understood as follows. The linear wire antenna, which induces the 
wave magnetic field in the rry-plane of variation, essentially carries current (J) along 
positive z-direction i.e., out of plane of the page. A constant magnetic field (Bo) is 
already present in the y-direction. This ambient magnetic field interacts with the 
antenna current, and gives rise to a net J x Bo force, which is directed along —x 
direction. Since the nonlincarities are at work in the large amplitude simulation, 
they exert the resultant Lorentz force (i.e., J x B force) on the wave field. The 
induced magnetic field of the wave thereby get influenced by the Lorentz force and 
consequently show a tendency to move in the direction of this force. Furthermore, a 
combined action of the wave motion along y-direction and the push in —x direction, 
overall results in slight distortion of the wave field topology. This is indeed what we 
observe in the simulation also (see Fig. (8) at t = 20). Our numerical investigations, 
thus demonstrate that nonlinear interactions significantly influence the propagation 
properties of whistler waves in EMHD. 

Nonlinear evolution of whistler waves is widely observed in many ionospheric 
and magnetospheric phenomena. For instance, whisthcr mode has been reported 
in the context of equatorial region of Earth's magnetosphere based on a modula- 
tional instability analysis [27], where it was shown that amplitude modulation in 
whistlers develops on much faster time scale than the wave transit period. Nonlinear 
interaction of whistlers is also likely to occur with magnetosonic fluctuations [28]. 
On the other hand, the nonlinear whistlers can potentially excite electrostatic low- 
frequency collisional gradient drift modes that are relevant to understanding the 
density irregularities in the lower part of the ionosphere [29] . Additionally, electro- 
magnetic whistler wave, propagating parallel to an external magnetic field, are also 
reported to decay into another circularly polarized electromagnetic wave and an ion 
acoustic wave in a homogeneous plasma [30]. This result remains unaltered under 
the combined effects of the relativistic-mass and ponderomotive-force nonlineari- 
ties [31]. Whistler waves in inhomogeneous media play a critical role in ionospheric 
density striations [32]. Our simulations of nonlinear evolution of whistler waves, 
relevant to understanding a number of features described in the above work, thus 
reveal that it is the hall force corresponding to the J x B term in electron mo- 
mentum equation that plays a crucial role and it is responsible for governing the 
nonlinear interactions associated with the electron whistler modes. 
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8. Conclusion 

In conclusion, we have carried out a detailed study of whistler wave propagation 
with the prime objective to understand their interaction and the nonlinear prop- 
agation characteristics. For this purpose a variety of initial configuration were en- 
visaged and their evolution studied. Whistler waves being dispersive, in the linear 
small amplitude limit different modes of the wave packet propagate with different 
phase velocities. However, the entire wave front moves with the typical group veloc- 
ity. The shape of the wave front is finally governed by the group speed. The linear 
propagation characteristics are found to be in agreement with the experimental 
observations. In the nonlinear regime we show that the propagation is significantly 
altered. This indicates the importance as well as the subtle role of nonlinearity 
in the context of whistler wave propagation thereby countering earlier speculations 
and statements about robustness of these modes against nonlinear effects. Our stud- 
ies presented in this paper can be useful in understanding a number of processes 
associated with the linear whistler waves, their mutual interaction and nonlinear 
feature in the context of interplanetary ionospheres and radiation belts. 
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